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ABSTRACT 



Circumbinary (CB) disks have been proposed as a mechanism to extract or- 
bital angular momentum from cataclysmic variables (CVs) during their evolution. 
As proposed by Taam & Spruit, these disks extend outwards to several a.u. and 
should be detected observationally via their infrared flux or by absorption lines 
in the ultraviolet spectra of the CV. We have made use of archival HST/STIS 
spectra as well as our own near-IR imaging to search for observational evidence 
of such CB disks in seven CVs. Based on the null result, we place an upper limit 
on the column density of the disk of Ah ~ 10 17 cm -2 . 
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1. INTRODUCTION 

Cataclysmic variables (CVs) are mass transferring binary star systems consisting of 
a white dwarf primary star and red dwarf secondary star. The evolutionary path from 
detached binary star system to close, interacting binary star system is greatly affected by 
orbital angular momentum loss (OAML), of which there must be significant amounts in order 
to transform the originally long period progenitor into an interacting binary (see Howell et 
al. 2001). A large amount of the OAML occurs during the common envelope (CE) phase of 
the CV evolution. 

Although the exact mechanism for OAML is still a matter of debate, several theories 
have been advanced to explain the observed orbital period-mass transfer distribution of CVs. 
Generally, magnetic braking is invoked for long period systems (P or b > 3 hr), and although 
these models are marginally successful in explaining the observed distributions, they predict a 
strong correlation between mass transfer rate (M) and orbital period, along with an excessive 
number of ultra-short period, very low luminosity systems; neither the stringent M — P or b 
correlation, nor the preponderance of low luminosity systems are observed. The discord 
between theory and observation could be due to a number of different factors, the most 
intriguing of which is that magnetic braking may not be solely responsible for OAML, but 
rather angular momentum may be extracted from the binary system via a circumbinary 
(CB) disk (Spruit & Taam 2001; Taam & Spruit 2001; Dubus et al. 2002; Taam, Sandquist, 
& Dubus 2003). The gravitational torque exerted by the CB disk on the binary will cause 
significant OAML leading to an accelerated period evolution and allowing the secondary star 
to completely dissolve within a finite time (less than a Hubble time). These two consequences 
are the attraction of the CB disk theory, as they would cause a broad range of M for any 
given period and remove the need for a large population of faint CVs. 

If CB disks exist, they likely begin formation during the CE phase and acquire material 
via mass outflow from the white dwarf (novae) or accretion disk (winds and outbursts) (Spruit 
& Taam 2001). For the CB disk to become a significant source of OAML, the mass input rate 
into the CB disk must be 10~ 5 to 0.01 of the material transferred from the secondary to the 
white dwarf throughout the CB disk lifetime. CB disks are assumed to be geometrically thin 
with a pressure scale height to radius ratio of H/R ~ 0.03, have a viscosity a < 0.01, and 
of Keplerian structure. Dubus et al. (2002) have proposed that the inner edge is truncated 
due to tidal forces from the binary, has a radius near 1 x 10 11 cm, effective temperature 
T ~ 3000 K, and is optically thick, while the majority of the outer disk is cool, T < 1000 
K, optically thin, and can extend out to several a.u. Continuum contribution comes from 
dust emission, with opacities due to dust grains, and includes contributions from iron-poor 
opacities at T < 1500 K where dust particles exist. Based on models of protoplanetary disks 
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(e.g., Bell et al. 1997), Dubus et al. (2002) find that their CB disks reach a maximum vertical 
height above the plane of the disk at a radius less than the outer edge of the disk. However, 
the exact physical height of the outer edge of the CB disk is not known (R. Taam, private 
communication). Figure 1 shows a schematic model of a CB disk based on specifications 
given in Dubus et al. (2002). 

Two properties of the proposed CB disks will manifest observationally; the large extent 
of the CB disk and its surface mass density, £ ~ 10 2 — 10 4 gcm~ 2 , should make the disk a 
source of "ISM-like" absorption lines appearing at the systemic velocity of the system, and, 
if the CB disk emission (from the warm inner edge) can be considered to be a sum of black- 
bodies, the flux from the CB disk will peak near 3 /j,m (Dubus et al. 2002). These properties 
lend themselves to detection in two different wavelength regions. Flux measurements via 
near-IR (1 — 10 /iin) imaging of face-on systems (offering a full view of the CB disk) may 
detect the IR flux of the CB disk. A second detection technique is to search for the ISM-like 
absorption lines produced by the CB disk in UV spectra of the binary systems. Hubble Space 
Telescope/ Space Telescope Imaging Spectrograph (HST/STIS) data are ideal for this as the 
high spectral resolution will enable detection and measurement of ISM-like lines whose radial 
velocities may then be compared to the space motion of the binary. This technique is best 
done for CVs with inclinations that match the opening angle of the CB disk, as the binary is 
then observed through the plane of the CB disk. The theoretical work (e.g., Taam & Spruit 
2001) predicts that the higher mass transfer rate CVs are the best candidates for a CB disk, 
as material may continuously be entering the CB disk, and may still retain some vertical 
extent above the disk (i.e., the CB disk has not had time to "flatten out"). Thus, we have 
begun an observational search for CB disks with long period nova-likes and old novae. 



2. FEASIBILITY OF CB DISK DETECTION 

While we have chosen CV candidates that should be ideal for detecting the proposed CB 
disks, the feasibility of disk detection depends quite a bit on the unknown CB disk structure. 
Here we discuss four different models for the disk structure - isothermal, isentropic, constant 
density scale height, and constant pressure scale height - and their significance for detecting 
a CB disk via their absorption signatures. 
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2.1. Key Assumptions 



The CB disk is expected to extend for several orders of magnitude in barycentric dis- 
tance, R, from the binary (Spruit & Taam 2001). The dependence of the vertically-integrated 
surface mass density, S, on radius, R, is taken from Taam & Spruit (2001), particularly Fig- 
ure 8 in that paper, from which we approximate 



Here Si n = 10 4 gcm~ 2 , and R extends from i?i n to R ut, which we take to be 10 11 cm 
and 10 14 cm, respectively (Dubus et al. 2002). Unlike a Shakura-Sunyaev disk, the opening 
angle, H/R, is expected to be roughly constant. We choose H/R = 9 opcn ~ 0.03 ~ 1.7°, 
taken from Figure 2 of Dubus et al. (2002). 

This value for ^open, as well as the relation for E in eq. (1), are assumed to be the same 
for each CV system. The difference between each system is then contained wholly in the 
inclination, i, or equivalently, the parameter # v icw, the viewing angle, where 90° — i — 6> vicw . 
Note that the smallest viewing angle in our sample occurs with EX Hya, where # v i ew = 
12° ± 1°, which corresponds to a height above the midplane of roughly z/H ~ 7, where 
z = R sin ^view — -Review is the height above the disk. 



In order to estimate the density above the midplane, we need a model for the vertical 
density profile of the disk. If we are within one pressure scale height of the midplane, then 
the uncertainty in the density is relatively small, but the theoretical uncertainty becomes 
much greater for higher \z\. Various processes, such as turbulent transport, transport by 
sound waves, and so forth, affect the vertical entropy profile, which in turn determines the 
vertical density profile. In hot disks, there is the further potential complication of winds 
being driven off the disk or a hot corona. Even in relatively cool disks there may also be 
photoevaporation off the disk by the illumination from the central binary. We leave these 
considerations to future theoretical work on the subject of CB disks. 

In calculating a vertical density profile for a disk, it is standard to assume a plane-parallel 
atmosphere. Note that in this case there is a vertical gravity in the disk that increases linearly 
with z: 




(1) 



2.2. 



Vertical Disk Profile 




(2) 
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where Qk is the Keplerian angular velocity, so that the equation of vertical hydrostatic 
equilibrium is 



This equation can not be solved without some additional physics to inform us of the vertical 
entropy profile, such as the local heating rate, q+, due to turbulent dissipation, and the 
radiative and turbulent transport of energy out of the disk. Some insight into this problem 
is provided by theoretical work (e.g., Cassen 1993), but uncertainties remain, especially at 
altitudes of several scale heights. 

First let us comment on the disk thickness (or half-thickness), H. From a disk theorist's 
perspective, H is usually taken to be a rough quantity that comes about in basic mixing- 
length-theory arguments. Given a location in the disk, H is determined by the sound speed, 
which in turn is determined by the temperature. All of this is tied together by standard 
thin-disk theory, in which quantities are taken to be functions of R only, so that the vertical 
structure appears only through some rather rough approximations. 

It is common to refer to this H as "the" pressure scale height, but when considering the 
vertical structure in greater detail, this definition is no longer sufficient. The local vertical 
pressure scale height, hp, defined as hp 1 = —(d z P)/P, may vary significantly with height, 
z, above the midplane; in fact it becomes formally infinite at the midplane. So, from the 
operational standpoint of an observer, there may be some significant ambiguity in defining 
H, which we need to address. While there have been theoretical studies of the vertical 
density profile for certain disks, we are not aware of any theoretical treatment of the vertical 
density profile for CB disks that extends to several scale heights above the midplane. What 
we offer here is a few simple models for the disk structure. Hopefully future theoretical and 
observational work will further constrain these models. 

Mathematically, one option for defining H is that it is the normalized second moment 
of the density, 



An alternative definition of H that is perhaps of more use to the observer is in terms of 
the photospheric surface r = 2/3. We also note that Bell et al. (1997) define H in terms of 
the "midplane pressure scale height", c s /Q. For the isothermal disk and the isentropic disk, 
we take H = Hi as defined above, and for the constant density scale height and (modified) 
constant pressure scale height models, we take H to be equal to the local density and pressure 
scale heights, respectively. 



d z P = -pQ 2 K z. 



(3) 




(4) 
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2.3. Model A: Isothermal Disk Atmosphere 

For a disk that is vertically isothermal, at a given radius, R, we have p ~ P. The 
equation of hydrostatic equilibrium, eq. (2), then gives d z p ~ —pz, which has solution 
p = po ex P ( — az 2 ). The local pressure scale height is hp 1 = 2az. If we set hp = z = H , we 
find 2aH 2 = 1, so that 

p = p exp(-z 2 /2H 2 ). (5) 

It is easy to verify that the H in this solution coincides with the definition for H\ given 
above, so setting hp = z = H is appropriate. Note that the local pressure scale height h P 
(which in this case is also the local density scale height) is then 

tt2 

hP = fj-, (6) 

\z\ 

i.e., it diverges at the midplane, is equal to H at z — ±if, and asymptotes to zero at infinity. 
The corresponding surface density, S, for this solution is 

r+oo 

S = / p(z) dz = V2^p H. (7) 



In terms of E and the height above the disk, z, the density becomes 

S 



HV2tt 



exp(-z 2 /2H 2 ). (8) 



2.4. Model B: Isentropic Disk Atmosphere 

For an isentropic disk, we begin with the equation of state P = Ap J , where A = P / p$, 
and Pq and p are fiducial values of the pressure and density, respectively, taken at the 
midplane. From equation (3) we have that 



which may be integrated to find 



dp Q| 



z dz. 



(9) 



P = Po 



7-1 



(10) 



A very notable aspect of this solution is that p = for \z\ > z max for some z max . Note 
that £ max may be written as a function of A, Qk, and 7, although we do not reproduce the 
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solution here. Assuming that 7 = 5/3, we find that 



(11) 



We wish to know the H that corresponds to the above solution. Here the matter is not 
so simple as for the isothermal case, and we must take care to define H appropriately, as 
discussed above. In the case of H defined in terms of the normalized second moment of the 
density, i.e., Hi, we find 

Hi = ~ (2.50)-^ max , (12) 
V6 

so that, if we set H = Hi, we may now write 

£ 8 



' H3nVG 
for \z\ < H\/Q, and p = otherwise. 



1 - 



hVq 



3/2 



(13) 



2.5. Model C: Constant Density Scale Height Disk Atmosphere 

The constant density scale height disk is rather trivial. If we define H as the local 
density scale height, 

K 1 = -^ = H, (14) 
P 

we find that 

P = ^exp(-\z\/H). (15) 



2.6. Model D: Constant Pressure Scale Height Atmosphere with Isothermal 

Core 

The constant pressure scale height atmosphere would seem to be an attractive model. 
However, it presents an unphysical solution. Namely, if we set 

h? = ~?f = H, (16) 



we find through equation (3) that 



p~^exp(-H/#), (17) 



\z\ 
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so that that density diverges at the midplane, and moreover the vertical integral of the 
density, i.e., E, diverges as well. 

As an alternative, we consider here an isothermal vertical profile near the midplane (the 
"core"), with a constant pressure scale height solution "glued" onto the isothermal solution 
at \z\ = H. This has the somewhat unsatisfactory result that the vertical derivative of the 
density is discontinuous at \z\ — H, i.e., there is a "kink" in the solution, but this way the 
divergence mentioned above is avoided. We then have 

P\z\<h = p exp(-z 2 /2H 2 ) (18) 

and 

P\ Z \ >H = Poe 1/2 ^e*v{-\z\/H) (19) 

where p is determined by the error function and the exponential integral, yielding the 
numerical value 

p ^ (0.4107)^. (20) 



We may then write 



and 



p\ z \ <H = (0.4107)-| exp (-z 2 /2H 2 ) (21) 
P\ z \>h = (°- 6773 )ff| exp(-|z|/if). (22) 



2.7. Column Densities for the Models 

It is important to note that, in all cases, the solution p(E, H, z) may be written in the 

form 

p(X,H,z) = ^f(\z\/H) (23) 

for some function /(). 

The column number density of hydrogen, Nu, is 

N u = J n H d£ (24) 

where d£ is an infinitesimal displacement along the line of sight to the binary. We assume 
here that the hydrogen mass fraction of the disk is near unity, so that n H — p/ra^. 

In peforming the integral, the key assumptions are: a constant flare angle H/R with 
radius R, the expression given above in eq. (1) for the dependence of E upon radius, and the 
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functional form in eq. (23) for the density. Furthermore, we assume that the density drops 
immediately to zero when R is outside the range (i?i n , -Rout)- Note that d£ = sec 9 vicw dR ~ 
dR. We therefore have 

i rRout 

N u = / pCZ, R, z) dR (25) 

m H J Rin 

The argument of /() that appears in the expression for p in eq. (23), namely \z\/H, is 
a constant, equal to tan# vicw / tan# opcn v i CW /9 opcn . The choice of E ~ R^ 1 is rather 
fortuituous from the point of view of performing this integral. Given that R ou t Rm and 
using the small-angle approximation, we find 

N H = — ^E in . (26) 

n^Yl "open 



where 77 = 9 vicw /9, 



open • 



The numerical value of the column density for our four models may then be written as 
follows. For Model A, the isothermal model, 



27 . / ^in \ ( ^open ^ f f] 

10 4 gcm" 2 J \T7° 



N H = (8.03 x 10^)1 ~ m _ 2 ^ «P -7T . (27) 



For Model B, the isentropic model, we find 



J V H = (6.98X10-)^ I5I |^)(^) (l-i.) . (28) 
For Model C, the constant density scale height atmosphere, 



-1 / o\ 3/2 



^ M *^\W^)\$F) eXpH "' ) ' <29> 

and for the constant pressure scale height atmosphere with isothermal core, 

N H = (8.27 x 10 27 ) ( ( ^ " exp ( -t) for \rj\ < 1 (30) 



10 4 gcm- 2 J V 1-7° 

N H = (1.36 x 10») (jgj^zj) (jif) axp(-N) for M > 1. (31) 

Figure 2 displays these values as functions of z/H(— 8 vicw /9 opcn ). 
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2.8. Discussion of the Solutions 

As may be seen in Figure 2, the expected density profile for small values of \z\/H 
{\z\ < 2H) does not depend too much upon the model chosen. Essentially, p is expected to 
be of the order of £/if , dropping gradually as \z\ increases, and this is what is seen. 

The situation changes quite dramatically, however, when we consider much larger \z\, 
as can be seen in Figure 2. As discussed above, the density in the adiabatic solution drops 
to zero at roughly \z\/H ~ 2.50, for 7 = 5/3 and H = Hi. The density remains rather high 
in the two constant scale height solutions, dropping exponentially or close to exponentially 
depending upon whether it is the density scale height or the pressure scale height that is kept 
constant. The isothermal solution is a sort of compromise between these extremes, dropping 
as a Gaussian. It is unlikely that a CB disk with a structure described by Model B would be 
detectable in any CV system that does not have % ~ 90°. CB disks with structures described 
by models A, C, and D, however, would be detectable in systems with lower inclinations. 

3. OBSERVATIONS 

As this project is an initial investigation for CB disks around CVs, we chose to look at a 
concise sample of objects that met our criteria: bright, high mass transfer systems (including 
the old nova, DI Lac); known inclinations and systemic velocities; and existing HST/STIS 
spectra that would be retrievable via the Multimission Archive at Space Telescope (MAST). 
The six objects we chose for spectroscopic investigation are listed in Table 1 and the data 
sets utilized are given in Table 2. Each HST/STIS observation employed the E140M grating 
in Echelle mode, which spans the wavelength range 1150 — 1735A, and has a resolving power 
of R = 45, 800, or ~ 0.02 A pixel -1 throughout each spectrum. The data were retrieved from 
the MAST archive and processed using IRAF/STSDAS software. We used IRAF/SPLOT 
to analyze the spectra. 

Four of the spectra are of excellent quality: IX Vel and V3885 Sgr both have a signal- 
to-noise ratio of ~ 30 in the continuum, and QU Car and EX Hya have S/N ~ 20 in the 
continua. The spectra of DI Lac and TZ Per are slighty noisier, with signal-to-noise ratios 
of 5 and 10, respectively. The DI Lac spectrum is too noisy shortward of 1200A to properly 
identify any spectral features below this wavelength. 

We also present near-IR JHKL' band photometric data of the nova-like CV, V592 Cas- 
siopeia. V592 Cas was observed in JHK at the Wyoming Infrared Observatory (WIRO) 
in 1997 (Ciardi et al. 1998) and in L' at the Fred Lawrence Whipple Observatory 1.2 m 
telescope on Mt. Hopkins for a total of 5 minutes on 1999 January 1 UT. The JHKL' data 
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were obtained as part of a campaign to study IR properties of CV secondaries. Additional 
JHK observations of V592 Cas have been secured at the United Kingdom IR Telescope on 
Mauna Kea and are in agreement with the WIRO data. 

4. ANALYSIS AND DISCUSSION 

For each spectrum, we identified absorption lines of ionic species commonly associated 
with the ISM and CV systems, such as C, N, Mg, Si, and S. We considered absorption lines 
that are above the 2a error level in flux and that have full width at half maximum (FWHM) 
values < 0.2 A, as absorption lines with larger FWHM values would not be associated with 
the slowly rotating CB disk or the ISM. Absorption lines were fit with a single Gaussian 
function to determine line centers and equivalent widths (EWs). Some absorption lines are 
composed of well-separated multiple absorption components (for example, the absorption 
lines of QU Car); for these we measured the line centers of the components individually. 
At the resolution of the E140M grating, our spectra are not sufficient to resolve closely 
spaced ISM features; for absorption lines containing unresolvable components, we measured 
the lines with a single Gaussian. Radial velocities for all measured absorption components 
were then calculated and are given in Table 3, along with the EWs of the absorption lines. 
The measurement errors on the radial velocities for narrow lines that are unblended and 
unsaturated are within the velocity resolution of the spectrum, which is typically ~ 3 km s -1 . 
Lines with EW > 100 mA are likely saturated or close blends of two strong components. In 
such circumstances, it is difficult to assign a meaningful centroid velocity to these features; 
we therefore assign uncertainties of ±7kms _1 to the radial velocities calculated for these 
lines. Absorption lines of excited state ions are noted in the table by an asterisk. 

4.1. ISM Lines of Sight 

Many of the identified absorption lines will be due to absorption by the ISM, in particular 
the local interstellar cloud (LIC); the warm, partially ionized cloud in which the Sun resides. 
The heliocentric velocity of the LIC is 25.7kms _1 towards the Galactic coordinates / = 186?1, 
b = — 16? 4 (e.g., Wood et al. 2002). Evidence of the LIC is frequently observed over a large 
fraction of the sky as a distinct range of velocities associated with ISM absorption components 
present in high dispersion UV spectra of nearby stars (e.g., Lallement & Bertin 1992; Holberg, 
Barstow, & Sion 1998; Wood et al. 2002; Kimura, Mann, & Jessberger 2003). In addition 
to the LIC, many other velocity components are also observed in certain directions. The 
LIC and these other components are commonly discussed in terms of an assemblage of warm 
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interstellar clouds within 30 pc of the sun. A good up-to-date and lucid description of these 
ISM velocity components is provided in Frisch, Grodnicki, & Welty (2002, hereafter FGW). 
In Table 1, we give the velocity of the LIC in the direction of each of our candidate objects. 

The LIC velocity is easily recognized in the histograms we have created for the ISM-like 
absorption lines found in the spectra of the six CV systems we studied. These histograms 
are shown in Figures 3-8. For four of the six objects, DI Lac, IX Vel, TZ Per, and QU Car, 
the LIC velocity is at or near the peak velocity in the histograms. EX Hya and V3885 Sgr 
do not follow suit; the peak in the EX Hya histogram is ~ Tkms" 1 from the LIC velocity, 
and neither of the peaks in the V3885 Sgr histogram are aligned with the LIC velocity. 

The line of sight to DI Lac lies in a direction expected to show LIC components. How- 
ever, the projected LIC velocity very nearly coincides with the systemic velocity of the 
system, making it difficult to observationally distinguish either (see Figure 3). Complicat- 
ing matters is the presence of lines due to interstellar C I, an ion not present in the LIC 
because it is easily destroyed by the ambient stellar UV radiation field. The presence of C I 
lines indicates that the line of sight to DI Lac passes through a diffuse cloud with sufficient 
density to preserve the C I ion. The presence of an unusually dense component of the ISM 
is also indicated by the very strong and saturated lines of S II. The mean velocity of the C I 
lines and the S II lines (— l.TSkms" 1 ) contributes significantly to the peak in the velocity 
distribution adjacent to the predicted LIC velocity. It would be tempting to associate these 
lines with a CB disk, however, the velocities are not consistent and the physical conditions 
in the CB disk are not appropriate to the existence of the C I ion. Rather, it appears that 
the line of sight to DI Lac passes through a relatively dense interstellar cloud, unrelated to 
the CV system. There is no evidence in the spectrum of DI Lac of lines arising from a CB 
disk. 

The line of sight to EX Hya lies in a direction that shows no LIC components. The 
peak of the radial velocity histogram (Figure 4) does, however, coincide with the expected 
velocity of the Filament structure described in FGW, consistent with several stars that lie 
in the same general direction as EX Hya. There is no evidence for a CB disk and the line of 
sight to this star is consistent with the ISM only. 

For IX Vel and QU Car, the line of sight to each object lies in a direction expected 
to show LIC components and the peaks of the histograms coincide with the expected LIC 
velocity (see Figures 5 and 6). The line of sight to IX Vel is consistent with the ISM only. 
The spectrum of QU Car shows several relatively weak C I lines, indicating the line of sight 
passes through a relatively dense interstellar cloud. There is no indication of CB disk related 
absorption for QU Car. 
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The line of sight to TZ Per lies in a direction expected to show LIC components, however, 
the histogram of velocities (see Figure 7) shows a broad distribution of velocities. Like the 
spectrum of DI Lac, the spectrum of TZ Per contains C I lines as well as relatively strong S II 
lines, consequently a similar interpretation - a line of sight through a dense interstellar cloud 
- holds for TZ Per. Many of the features in the TZ Per spectrum are either saturated or quite 
strong, indicating that several unresolved velocity components are present and contributing 
to the broad distribution of velocities. There is one C I line appearing near the systemic 
velocity, however, this is not sufficient evidence to support absorption from a CB disk and 
as discussed above, the existence of this ion in a CB disk is not probable. 

The line of sight to V3885 Sgr contains two distinct, well separated velocity components 
(see Figure 8). The component at +2.5 km s" 1 can be confidently identified with the South 
Polar Cloud (see FGW). The second peak near 20kms _1 corresponds to no previously de- 
scribed interstellar cloud in the local ISM. There is no clear contribution from the LIC or a 
CB disk. The line of sight to this star is consistent with the ISM only. 

4.2. Absorption Line Velocities 

Figures 9-14 display line velocity profiles for select absorption lines from the spectra 
of each of the six CVs. On each plot, the vertical dashed line represents the LIC veloc- 
ity towards the object, and the vertical dotted line represents the systemic velocity of the 
system. For systems that have systemic velocities with large errors, the dotted line repre- 
sents the measured value. We selected a sampling of absorption lines to illustrate the line 
structure and velocity behavior for each system. Shown are the common absorption lines 
of the N I AA1199.55, 1200.22, 1200.71 triplet, lines from several C I multiplets, the S II 
AA1250.58, 1253.81, 1259.52 triplet, lines from several Si II multiplets, along with the Fe II 
A1608.46 and Al II A1670.79 absorption lines (excited transitions are noted with an asterisk). 

Figure 9 displays absorption lines found in the spectrum of DI Lac. Each line exhibits 
a single velocity structure, however, because the LIC velocity and systemic velocity are 
separated by less than the velocity resolution of the spectrum, it is not possible to identify 
the absorption lines with either velocity. As discussed in the previous section and based on 
the results from the remaining objects of our study (see below) it is likely that the absorption 
lines in the DI Lac spectrum are ISM in nature. The 1910 nova eruption of DI Lac implies 
the presence of nova shells, which may be observed via absorption lines, however, these lines 
should be blue-shifted with respect to the systemic velocity. 

The absorption lines found in the EX Hya spectrum (shown in Figure 10) have a single 
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velocity structure. They are shifted ~ 7 — lOkms -1 with respect to the LIC velocity of 
— 14.8 km s -1 , but are aligned with the Filament structure discussed in FGW. The line ve- 
locities of IX Vel (Figure 11) and QU Car (Figure 12) are aligned with the LIC velocity. In 
each case, the systemic velocity of the system is clearly not associated with the absorption 
line velocities. 

TZ Per (Figure 13) is another system for which the LIC and systemic velocities fall 
within the velocity width of the absorption line profiles. The two velocities are separated by 
more than the velocity resolution of the spectrum, therefore the absorption line velocities 
may be identified with the LIC velocity, as it is more closely aligned with the absorption line 
centers (the LIC velocity alignment is also seen in Figure 7). 

The absorption lines in the V3885 Sgr spectrum (Figure 14) have two velocity compo- 
nents; one that is aligned with the South Polar Cloud (FGW) and one that does not match 
any known ISM velocity system. Sagittarius is in the direction of the galactic center, so it 
is likely that different clouds along this line of sight are causing absorption along the line of 
sight to V3885 Sgr. 

4.3. Infrared Observations 

Although the outer regions of the proposed CB disks are cool and optically thin and not 
a source of IR blackbody flux, the warm, optically thick inner regions of the disk will make 
a contribution to the IR flux of a CV system. Given that the inner region of the CB disk is 
predicted to have a temperature of up to ~ 3000 K and is optically thick, we approximated 
a mean temperature of the CB disk that agrees well with Figure 2 of Dubus et al. (2002) 
to create a blackbody model representative of the CB disk. Determining this temperature 
is important, as the peak of the blackbody curve of the CB disk should cause an observable 
'bump' in the spectral energy distribution (SED) of the CV. 

V592 Cas has an inclination of i = 27°, and is therefore ideal for detecting the full 
CB disk via IR observations. In Figure 15, we plot JHKL' points for V592 Cas, along with 
blackbody curves representing the secondary star for the system and a steady state accretion 
disk. Also plotted is the sum of the IR flux from all components of the CV system, including 
the flux contribution expected from a CB disk; it is a sum of blackbodies and peaks near 
3 fim. This model summed spectrum overestimates the IR flux observed for V592 Cas; the 
data points do not fall along the model curve that includes the flux contribution of the CB 
disk, but rather follow the blackbody curve of the steady state accretion disk. It is apparent 
that flux from a CB disk is not evident in the IR data of V592 Cas. 
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4.4. Observations of the ISM 

Most of the spectral features observed in the six CV systems studied can be reasonably 
assigned to the LIC or other recognized velocity structures in the local ISM near the Sun. 
However, in several CVs (DI Lac, QU Car, and TZ Per), we do observe unusually strong 
lines due to S II and features due to C I, indicating that the line of sight to these stars 
passes through a dense, warm ISM component where C I is shielded from ionizing radiation. 
It is regarded as unlikely that these regions are associated with either CB disks or other 
gas related to the CV, as the observed velocities of these components are distinctly different 
from the CV system velocity. 

4.5. Evidence for CB Disks? 

In the spectroscopic and photometric data presented here, we have found no convincing 
evidence to support the theory of CB disks around CVs. The ISM-like lines appearing in the 
spectra of the six CVs have velocities consistent with the ISM in the direction of the object. 
The two cases in which the systemic velocity fell within the velocity width of the measured 
absorption lines (DI Lac and TZ Per) were complicated by the fact that the systemic velocity 
was close to the LIC velocity and that the identified lines were unlikely to occur in a CB 
disk. Based on the results of the other four objects, the match of the systemic velocity is 
due purely to its proximity to the LIC velocity. Evidence for a CB disk was not found in the 
JHKL' photometric data of V592 Cas. There is no excess of flux around the expected peak 
wavelength of 3 fxm, and models of a CB disk overestimate the flux measured in V592 Cas. 

In the absence of any absorption features clearly associated with CB disks, the upper 
limits for the disk material may be derived based on 3a upper limits for the EWs of N I, 
Si II, O I, and S II, evaluated at the systemic velocity in the spectrum of each CV (except for 
DI Lac, which exhibits some lines near the systemic velocity). The 3<r EW values for all lines 
measured are very low - they range between 5 mA - 25 mA - so we can use the linear portion 
of the curve of growth to determine the species column densities: Ni = 1.129 x I0 20 (^ i )(-rr)- 
Wavelengths and / values are taken from Morton (1991). To calculate equivalent H column 
densities, we use typical abundances of each species observed for the warm phase of the ISM, 
Ai, such that Nu = jjS The best upper limits we calculate are from the N I and O I lines. 
For EX Hya we determine a value of JV H ~ 3 x 10 17 cm- 2 , for IX Vel N H ~ 1 x 10 17 cm~ 2 , 
for QU Car N H ~ 2 x lfFcm" 2 , for TZ Per N H ~ 2 x lfFcm" 2 , and for V3885 Sgr 
JV H ~ 2 x 10 17 cm- 2 . 

We have indicated the observational upper limit of Nr ~ 10 17 cm~ 2 as the dashed 
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horizontal line in Figure 2 for comparison of the observational numerical values with the 
theoretical disk models. The vertical dotted line at z/H ~ 7 represents the value of z/H for 
EX Hya, while the vertical dotted line at z/H ~ 17 represents the value of z/H for IX Vel, 
QU Car, and V3885 Sgr (the values for DI Lac and V592 Cas are past the limits of the plot). 
If CB disks exist, we can eliminate Models C and D as the proposed CB disk structure 
for the four CVs, as the model column densities are significantly above the detection limit 
and absorption by such disk structures would be detectable at the given inclinations. A CB 
disk with the structure of Model A would be detectable for the inclination of EX Hya and 
can therefore be ruled out, as no absorption from a CB disk is seen in EX Hya. Model B 
would not be observed in any of our CV systems, but we note that this model is the most 
conservative disk assumption of the four and we do not expect it to be physically realistic 
(i.e., we would expect the CB disk to have some form of an atmosphere). 

As this investigation has only considered a small sample of CV systems that could 
contain a CB disk, we plan to continue our photometric and spectroscopic investigation for 
CB disks. In particular, we plan to analyze additional archival STIS data and to continue 
to look for their signatures in the IR SEDs of candidate CVs. 

We would like to thank D. Ciardi for the use of his V band image of V592 Cas prior 
to publication. KEB would like to thank R. Taam for helpful discussions regarding CB 
disks. We also thank the anonymous referee for comments that greatly helped to improve 
this manuscript. NS acknowledges support from NSF/REU grant AST-9810770. SBH ac- 
knowledges support from NSF grant AST-9810770. 
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Fig. 1. — A schematic of a CV orbital plane showing the inner edge of a model CB disk. The 
inner disk starts near 10 11 cm, is warm and optically thick. The outer edge is cold, optically 
thin, and extends out to several A.U. 
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Fig. 2. — Vertical density profiles above the CB disk midplane for disk structure models A, 
B, C, and D. The vertical dotted line at z/H = 7 represents the value for EX Hya, while 
the line at z/H = 17 represents the z/H value for IX Vel, QU Car, and V3885 Sgr. The 
horizontal dashed line at logiVn = 17 represents the 3cr upper limit for the column densities, 
based on measurements of the HST spectra (discussed in a later section). 
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Fig. 3. — Histogram of the radial velocities of ISM-like absorption lines measured in the 
spectrum of DI Lac. Labeled are the systemic (7) and calculated LIC velocities, neither of 
which match the peak velocity in the histogram. 
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Fig. 4. — Histogram of the radial velocities of ISM-like absorption lines found in the spectrum 
of EX Hya. The LIC velocity is labeled at — 15kms _1 , while the velocity of the Filament 
structure discussed in FGW at — 9kms _1 is also labeled. 



- 23 - 



1 I — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — r~ 

LIC 

8 - 



6 - 




2 - 



q I I I I U I I I I I I I I I I I I I I U I I i I i I I U I I I I 

-15 -10 -5 5 10 15 

Heliocentric Velocity (km s ) 

Fig. 5. — Histogram of the radial velocities of ISM-like absorption lines found in the spectrum 
of IX Vel. The peak of the histogram matches well with the LIC velocity of 5.6 km s -1 . 
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Fig. 6. — Histogram of the radial velocities of ISM-like absorption lines found in the spec- 
trum of QU Car. The peak velocity of the histogram corresponds to the LIC velocity of 
— 6.1 kms -1 . 



-25 - 




Heliocentric Velocity (km s ) 



Fig. 7. — Histogram of the radial velocities of ISM-like absorption lines found in the spectrum 
of TZ Per. There is a grouping of lines with velocities around 5 — 15 km s" 1 - the LIC velocity 
falls at the high end of this range. One absorption line has a velocity that matches the 
systemic velocity (labeled as 7) at — 2kms~ 1 . 
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Fig. 8. — Histogram of the radial velocities of ISM-like absorption lines found in the spectrum 
of V3885 Sgr. The two peaks in the histogram, at ~ 4 km s _1 and ~ 22kms~ 1 , represent 
the two velocity components of the absorption lines. The peak at ~ 4kms~ 1 matches well 
with the South Polar Cloud. 
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Fig. 9. — Line velocity profiles for absorption lines from four C I multiplets and the S II 
triplet observed in the DI Lac spectrum (S/N ~ 5). The dashed line at 3.5 km s -1 represents 
the LIC velocity, while the dotted line at lkms" 1 represents the systemic velocity. 




Fig. 10. — Line velocity profiles for the N I triplet, lines from several Si II multiplets, the 
Al IIA1670 line, and the S II triplet found in the EX Hya spectrum (S/N ~ 20). The dashed 
line at — 14. 8 km s" 1 represents the LIC velocity. 
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Fig. 11. — Line velocity profiles for the N I triplet, lines from several Si II multiplets, the 
Al IIA1670 line, and the S II triplet found in the IX Vel spectrum (S/N ~ 30). The dashed 
line at 5.6 km s -1 represents the LIC velocity, while the dotted line at 48 km s -1 represents the 
systemic velocity. The velocity centers of the absorption lines clearly match the calculated 
LIC velocity. 
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Fig. 12. — Line velocity profiles for the N I triplet, lines from several Si II and C I multiplets, 
and the S II triplet appearing in the QU Car spectrum (S/N ~ 20). The dashed line at 
— 6.1kms _1 represents the LIC velocity, while the dotted line at — 84kms _1 represents the 
systemic velocity. The velocity line centers are in good agreement with the calculated LIC 
velocity. 
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Fig. 13. — Line velocity profiles for the N I triplet, lines from several Si II and C I multiplets, 
and the S II triplet appearing in the TZ Per spectrum (S/N ~ 10). The dashed line at 
14.7 km s" 1 represents the LIC velocity, while the dotted line at — 2 kins" 1 represents the 
systemic velocity. Both the systemic and LIC velocities fall within the velocity widths of the 
absorption lines, but it is the LIC velocity that most closely matches the line centers. 
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Fig. 14. — Line velocity profiles for the N I triplet, lines from several Si II multiplets, the 
Fe II A1608 and Al IIA1670 lines, and the S II triplet appearing in the V3885 Sgr spectrum 
(S/N ~ 30). The dashed line at — 20.0 km s" 1 represents the LIC velocity, while the dashed 
line at — 45kms~ 1 represents the systemic velocity. Each absorption line exhibits a two 
velocity component structure. 
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Fig. 15. — A blackbody model of the IR flux of V592 Cas, including the contribution from 
a CB disk, created using a steady state CV accretion disk code. The model for the CB disk 
is a "sum of blackbodies" and the flux peaks near 3 /im. Also plotted are the JHKL' data 
points for V592 Cas, which do not match the CB disk blackbody curve. 



Table I. Object Log 





I, b 


Inclination (z/H) 


LIC Velocity 


Distance 


Systemic Velocity 




Object (type a ) 


(dcg) 


(deg) 


(kms- 1 ) 


(PC) 


(kms- 1 ) 


References 


DI Lac (N) 


103.107, -4.855 


18 (42) 


3.5 


2500 


1 


1,2 


EX Hya (IP) 


303.186, 33.622 


78 ± 1 (7) 


-14.8 


64.5 ± 1.2 


-180 


3,4 


IX Vcl (NL) 


264.929, -7.890 


60 ±5 (17.6) 


5.6 


95 ± 12 


48 


5,6 


QU Car (NL) 


293.513, -7.717 


60 (17.6) 


-6.1 


2000 


-84 ± 20 


7,8 


TZ Per (DN) 


133.576, -2.777 




14.7 


80 


-2 ±2 


9,10 


V3885 Sgr (DN) 


357.477, -27.762 


50 - 70 (17.6) 


-20.0 


110 ± 25 


-45 ±4 


5,11,12 


V592 Cas (NL) 


118.603, -6.910 


28 ± 11 (36) 


9.9 


60 


21 ± 14 


13 



References. — (1) Moyer et al. (2003); (2) Kraft (1964); (3) Hellier et al. (1987); (4) Beuermann ct al. (2003); (5) Beuermann 
& Thomas (1990); (6) Perryman et al. (1997); (7) Gilliland & Phillips (1982); (8) Drew ct al. (2003); (9) Echevarria ct al. 
(1999); (10) Berriman, Szkody, & Capps (1985); (11) Cowley ct al. (1977); (12) Haug & Drechsel (1985); (13) Huber ct al. 
(1998) 

a N=nova, IP=intcrmediate polar, NL=nova-like, DN=dwarf novae 



Table 2. HST Observations Log 



Object 


Data Set 


Observation Date 


Grating 


Exposure Time 
(s) 


S/N 


DI Lac 


O5B607010 


2000 April 19 


E140M 


2127 


5 


EX Hya 


O68301010 a 


2000 May 


E140M 


15,200 


20 


IX Vel 


O5BI01010 


2000 April 3 


E140M 


1750 


30 


QU Car 


O5BI08010 


2000 July 21 


E140M 


2600 


20 


TZ Per 


O5B610010 


2000 July 3 


E140M 


2151 


10 


V3885 Sgr 


O5BI04010 


2000 April 30 


E140M 


2480 


30 



a The EX Hya data sets 068301010, 068301020, 068301030, 068302010, 
068302020, and 068302030 were used for the EX Hya measurements. 



Table 3. Identified Absorption Lines 







DI Lac 


EX Hy 


a 


IX Vel 




QU C 


&r 


TZ Per 


V3885 S 


err 
6 1 


Ion A r e S t 


RV 


EW a 


RV 


EW 


RV 


EW 


RV 


EW 


RV 


EW 


RV 


EW 




(A) 


( kms -1 ) 


(mA) 


( kin s~ ) 


(mA) 


( km s — ) 


(mA) 


( km s — ) 


(mA) 


( kms -1 ) 


(mA) 


( km s ^ ) 


(ml) 


N I 


A1159.82 














-12.28 


43 










C I 


A1188.83 














-8.00 


12 










C I 


A1190.02 


















13.80 


45 






Si II 


A1190.42 






-10.08 


122 


3.81 


93 


-10.53 


147 


6.24 


125 


1.00,23.68 


62,84 


Si II 


A1193.29 






-8.80 


103 


3.78 


112 


-8.02 


173 


15.06 


84 


6.09,21.17 


98,116 


C I 


A1193.03 














-2.99 


7 


13.81 


91 






Mn II 


A1197.18 














-5.53 


27 










Mn II 


A1199.39 














-3.04 


16 










Mn II 


A1201.12 














-8.05 


16 










N I 


A1199.55 






-5.60 


115 


1.31 


138 


-10.45 


182 


13.82 


110 


3.69,21.20 


93,117 


N I 


A1200.22 






-5.90 


105 


2.88 


126 


-11.38 


156 


16.32 


89 


6.20,23.70 


74,95 


N I 


A1200.71 






-7.67 


78 


3.75 


100 


-10.55 


131 


18.44 


127 


6.21,23.70 


74, 105 


Si III 


A1206.50 






-3.13 


68 


11.10 


73 


-3.10 


116 






11.28 


76 


Mg II 


A1239.93 














-5.79 


26 


11.48 


36 






Mg II 


A1240.39 














-0.96 


19 










SII 


A1250.58 


-1.64 


125 


-9.01 


11 


7.04 


43 


-6.81 


91 


8.16 


94 


5.81,26.20 


21,24 


SII 


A1253.81 


-9.03 


233 


-8.03 


35 


10.11 


27 


-6.11 


105 


6.51 


99 


4.67, 26.20 


24,48 


SII 


A1259.52 


-3.49 


182 


-7.98 


38 


8.12 


44 


-5.67 


111 


11.79 


109 


2.38,23.82 


26,56 


Si II 


A1260.42 






-8.86 


8 


4.99 


136 


0.73 


269 


18.84 


228 


7.16,23.82 


115,126 


CI 


A1260.74 


















9.18 


60 






CI 


A1276.48 


















-3.19 


32 






CI 


A1277.25 


-1.49 


96 










-1.31 


33 


13.97 


21 






CI* 


A1277.55 


















6.93 


50 






CI 


A1280.14 


1.47 


56 










1.33 


6 


13.98 


81 






PII 


A1301.87 


7.86 


64 










-3.83 


20 


9.41 


46 






01 


A1302.17 


-12.87 


196 


-9.31 


136 


0.43 


163 


-8.44 


195 


2.96 


289 


12.38,51.54 


195,17 



Table 3 — Continued 







DI Lac 


EX Hy 


a 


IX Vol 




QU Car 


TZ Per 


V3885 S 


gr 


Ion Arest 


RV 


EW a 


RV 


EW 


RV 


EW 


RV 


EW 


RV 


EW 


RV 


EW 




(A) 


(kms -1 ) 


(mA) 


(kms -1 ) 


(mA) 


(kms -1 ) 


(mA) 


(kms -1 ) 


(mA) 


( kms -1 ) 


(mA) 


( kms -1 ) 


(mA) 


c;; TT 




Q P.A 
— O.04 


OIK 


-8.28 


79 


5.01 


81 


— O.40 


loo 


11. 1 z 


1 a n 
14U 


3.20,23.90 


52,86 


i\ l ii 


\ 1 Qi 7 oo 

AlOl ( .ZZ 


— o.oo 


7/1 










— o.ZD 


on 










C T 

\j 1 


AloZo.oo 


1 01 

— l.zl 


91 










a no 
— 4.Uz 


ZZ 


y.oo 








CI* 


A1329.09 


0.0U 


/o 














11.81 


DO 






CII 


A1334.53 


-6.20 


44 


-10.79 


151 


1.93 


170 


-6.82 


213 


11.82 


244 


12.20 


214 


C II* 


A1335.66 


4.95 


166 


0.36 


28 


5.20 


54 


-10.36 


129 


7.79 


118 


14.97,32.94 


25,33 


CII 


A1347.24 














-6.28 


27 


7.41 


92 






Ni II 


A1370.13 


5.56 


177 






















SI 


A1425.03 


















3.69 


21 






SI 


A1425.19 


















5.80 


18 






Si II 


A1526.71 


-3.60 


525 


-8.21 


129 


3.67 


111 


-6.49 


172 


5.12 


153 


0.65,24.23 


73, 106 


PII 


A1532.53 


















10.45 


82 






CI* 


A1560.70 


-2.13 


110 






















Fell 


A1608.46 


5.63 


260 


-7.46 


42 


2.40 


49 


-7.45 


131 


6.99 


125 


1.95,22.47 


23,24 


Al II 


A1670.79 






-8.02 


138 


-11.25 


61 


1.80 


161 


16.79 


145 


1.06,20.81 


45,51 



a incs with EW > 100 mA are likely saturated or blended components; we therefore assign an error of ±7 kms 1 to the radial velocities of 
these lines. See text for details. 



